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In this paper, we study the equation of state admissible for a flat FRW models filled with a bulk 
viscous fluid by using the Lie group method. It is found that the model admits scaling symmetries 
iff the bulk viscous parameter 7 = 1/2. In this case, it is found that the main quantities follow a 
power law solution and in particular the bulk viscous pressure 11 has the same order of magnitude 
as the energy density p, in such a way, that it is possible to formulate the equation of state 11 = >cp, 
where >c £ R~ (i.e. is a negative numerical constant). If we assume such relationship we find again 
that the model is scale invariant iff 7 = 1/2. We conclude that the model accepts a scaling symmetry 
iff 7 = 1/2 and that for this value of the viscous parameter. If = >cp, but the hypothesis 11 = xp 
does not imply 7 — 1/2, and that the model is scale invariant. 



I. INTRODUCTION 



As it is known in General Relativity (GR) there is a dualism between spacetime and matter. While the structure 
of the spacetime is governed by the field equations the physical properties of matter are introduced through the 
\^ ' energy-momentum tensor attending to diverse physical considerations. Some of these considerations come from fields 
of the physic where gravity does not play any role and are assumed in GR. As the Einstein equations together Bianchi 
■ identities form an undetermined system of equations, it is necessary to introduce equation of state, some of them ad 
hoc, in such a way that the resulting system of equations may be integrated. Such system satisfy certain symmetries 
that form a group, the group of symmetries of the equations. 
] Collins ([!]) and later on M. Szydlowsky ([2], we follow closely this work) have used the inverse way in order to 
. determine the admissible equations of state for a system of equations under the restriction that this system admits 
T-H ' a determined group of symmetries i.e. one could a priori assume a symmetry group of the Einstein equations and 
. out of it deduce the condition of integrability which has the form of the equation of state. In this way the group of 
' symmetries of the Einstein equations correctly select physical meaningful equation of state. 
^ , Since the bulk viscous theory is constructed assuming phcnomcnological (ad hoc) equations of state and therefore 
O^' the equations depend on certain undetermined numerical constants, we are interested in determining the exact form 
I ' of such equations of state imposing the condition that the field equations admit a concrete symmetry. Therefore, in 
bJO' this paper, we show that the field equations of a cosmological flat FRW model filled with a bulk viscous fluid together 
^ ' with a suitable equations of state admits a certain Lie group of symmetries or, vice versa, the invariancc of the field 
• ^ . equations with respect to a given symmetry group singles out the corresponding equation of state. 

' The paper is divided as follows. In section 2 we outline the general field equations of our model i.e. a flat FRW 
bulk viscous fluid and without the cosmological constant. In section 3 wc deduce the second order ode that governs 
the model. This differential equation has been deduced without any assumption. One we have outlined the basic ode, 
using the Lie group technique we study this equation finding that only admits one symmetry but if we impose that 
the model admits the scaling symmetry then we find that this is only possible if the viscous parameter 7=1/2 (where 
the viscosity ^ has been introduced into the field equations through the ad hoc law ^ = k-^p'^ , and in particular, we are 
interested in determining the possible value(s) of the parameter 7). Once we have established that the field equations 
are scale-invariant iff 7 = 1/2 we are interested in finding the relationship between 11 and p (i.e. we are interesting in 
determining a new equation of state relating 11 and p). For this purpose we try to integrate the resulting ode under 
the restriction 7 = 1/2, following the standard Lie procedure but unfortunately we had not been able of obtaining 
any explicit solution. Nevertheless we have obtained the invariant solution (a particular solution) that induces the 
scaling symmetry finding in this way a concrete power law solutions for the main quantities of the model. In this way 
we arrive to the conclusion that 11 and p has the same order of magnitude and therefore wc find that 11 = xp where 
>«• is a negative numerical constant. These results are not new, they have already been obtained by several author 
using different methods and they will be commented in this section. 

As we have been able to determine a concrete relationship between 11 and p (under the scale-invariant condition 
7 = 1/2) in section 4 we investigate if this condition implies 7 = 1/2. For this purpose, imdcr this hypothesis, we 
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obtain a second order ode that describes all the model and when studying it with the Lie group method we find 
that such equation only admits scaling symmetries iff 7 = 1/2. We study the resulting ode finding the same results 
than in the above section. Nevertheless, the assumption 11 = >cp, allows us to obtain a complete solution to the field 
equations, this possibility will be show in subsection 4.2. 

In section 5 we again study the model as well as some of the odes that have been arising in the paper through a 
pedestrian method, Dimensional Analysis. In this section we shall show how this method works in order to obtain 
the "same" results but in a trivial way. We end by summarizing some results. 



II. THE MODEL. 



For a flat Friedmann- Robertson- Walker (FRW) Universe with a line element 

ds^ = c^df - f{t) [dx^ + dy"^ + dz^) , (1) 

filled with a bulk viscous cosmological fluid the energy- momentum tensor is given by (see [3]) 

Tl: = {p + p + \l)u,u'' -{p + n)5l (2) 

where p is the energy density, p the thermodynamic pressure, 11 the bulk viscous pressure (stress) and Ui the four 
velocity satisfying the condition UiU^ = 1. The field equations yield: 

2H' + SH^ ^-K{p + n), (3) 

^ up, (4) 

p' + iapH = ^mn, (5) 

n' + 1-^ = -ipH - in (?,H - , (6) 



where 



M^^^ = l + ^, a^(c. + l), (7) 

and where we are assuming the following phcnomcnological (ad hoc) equation of state (laws) for p, ^, T and r (see 
[3]): 

p = u:p, i^k^p\ T = Dpp^, r^^p-' = k^p'>-\ (8) 

where < w < 1, and > 0, Dp > are dimensional constants, 7 > and f3 > Q are numerical constants. 
Eq. {p — Lop) is standard in cosmological models whereas the equation for r is a simple procedure to ensure that 
the speed of viscous pulses does not exceed the speed of light. These equations are introduced without sufficient 
thermodynamical motivation, but in absence of better alternatives we shall follow the practice adopting them in the 
hope that they will at least provide indication of the range of possibilities. For the temperature law T — Dpp^ which 
is the simplest law guaranteeing positive heat capacity. 

For a detailed deduction of this model see R. Maartens ([3]) 



III. THE GENERAL EQUATION. 



Without any assumption the field eq. (3-6) may be expressed by a single one 



H" - KqH-^ [H' f + KiHH' + K2H'H^-'^^ + K^^H^ + K^H'^-'^^ = 0, 



(9) 
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where 

UJ 



i^o = = (1 + /?) , P 



w + 1' 



/o\l-7 ol-7 

since [k^] = L7-iAf i-tT^t-i, [k] = LM"!, hence: [JCs] = 

Taking into account the value of the constants Ki this equations yields: 

H" - WH-^ [H'f + ZHH' + K2H'H'^-'^'< + ^{uj - I) + KaH^-'^'< = 0, (10) 

and if we decide to make the following assumption k-^ = k = 1 then eq. (10) yields: 

H" - ( ^f^i±l^ {^H'f + 3HH' + 3i-^F'i/2-27 + 9 _ 1) 3^ ^) ^4-27 ^ q ^^i) 
\uj + l J ^ ' ' 2 ^ ' ^ ' 

We go next to study this equation under the Lie Group technique (see for example [4], [5] and [li]). For simplicity 
we have rewrite it in the following form: 

H" - AR-^ [H'f + 'iHH' + CH'H^-^'^ + MH^ + EH^-'^'^ = 0. (12) 

The standard Lie procedure brings us to obtain the next system of pdes: 

^hhH + AiH = 0, (13) 

Ar? + QinH^ + VhhH^ - 2^tHH^ - Ai^hH + 2^7^^/^^"'^ - 0, (14) 

2mtH - (ttH + 3r]H + S^tH^ + 3M^hH^ + C^tH^-^"' + SE^hH^'^^ + (2 - 2j) CtjH^-^^ - 2ATjt = 0, (15) 

3M7]H + MH^ (2Ct - m) + (4 - 27) ErjH^-^'' + CijtH^~^'' + (2^ - m) EH^"^'' + rjuH-^ + im = 0, (16) 

this system admits the following symmetry 

e = i,?/ = o^Xi = at, (17) 

since Xi span an algebra Li the equation cannot be completely integrated by the Lie group method. 
But if we try to check if the system admits a scaling symmetry 

e = t,77 = -i/, (18) 

this is only possible iff 7 = 1/2. Therefore 

C = ai + 5,77 = -aF, (19) 

is a symmetry of the ode iff 7 = 1/2, and where a and h are numerical constants i.e. a, 6 £ M. Hence 

Xi = du X2^tdt-HdH [Xi,X2]=Xi, (20) 

which span a solvable Lie algebra L2 of the type ///. 

At the same result have arrived for example A.A.Coley et al (["]), who study this model from a dynamical system 
approach. To apply this method they rewrite the field equations in a dimensionless way in such a form that this is 
only possible iff 7 = 1/2 and the viscous pressure and the energy density has the same order of magnitude as we 
will see in the next section. At similar conclusions R. A. Daishev and W. Zimdahl have arrived in ([S]), where these 
authors study this model from the homothetic (similarity) point of view i.e. they study when the field equations 
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remain self-similar. As we will see below we have obtained the same results but using the Lie group method. Finally 
Belinchon et al ([()]) have obtained the same results using the renormalization group approach. 
The canonical variables and the reduced ode that induces the symmetry Xi arc: 

yi^) = ^ = H, (21) 

y' ^ (Mx^ + i5.T4-27) y3 + (3^. + Cx^~2^) - A^, (22) 

X 

which is an Abel ode (see [10]). Without any assumption it is very difhcult to find any explicit solution of this equation 
and therefore a solution of eq. (12) 

A. The case with 7 = 1/2, scale invariant solution. 

As we can see, the ode (12) admits a scaling symmetry iff 7 = 1/2 and in this case such ode is reduced to: 

H" - WH~^ (H'f + (3 + Vs) HH' + (^^(uj-l) + ^iLu + 1) j = 0, (23) 

or equivalently 

H" - AH^'^ {H'f + BHH' + CH^ = 0, (24) 
where obviously it admits the symmetries 

Xi=dt, X2=tdt-HdH, [Ai,A2]=Xi, (25) 

which form a L2 algebra etc.. 

Symmetry Xi brings us to the following ode through the reduction (canonical variables) 



y{^) = ^^H, (26) 



y'=(^lic^-l) + ^{u; + 1) j x^y^ + (3 + x/3) xy^ W^, (27) 
y' = Cx^y^ + Bxy^ - A^, (28) 

X 

which is an Abel ode. This ode admits the following symmetry 

X = xd^~ 2ydy, (29) 
which is a scaling symmetry and it induces the following change of variables, 

r = x^y, s(r)=ln(a;), =^ a; = e*^''^ y=-^^^ (30) 

which brings us to obtain the next ode in quadratures 

(31) 



r (Cr2 + Br + 2- Ay 

and which solution is: 

Inr 1 In (Cr^ ^ Br + 2 - A) ^ ^'^^^^^ ^ ( Vb'+Ic(A-2) ^ 
'^'^-~T—2^2 A^2 {A-2),/WTW{T^ ^ ' 
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and hence in the original variables {x,y): 

_ ln(x2y) lln{Cx^y'^ + Bx'^y + 2-A) ^ Cretan /i 

lux — \- — , h ('J'j} 

A- 2 2 A- 2 (A-2)^S2+4C(A-2) 

therefore we have obtained the next ode in the (H, t) variables: 



In(fi) ^lniC[^f + B[^)+2-Aj Barctan/. 



lni? = ^^^T. ^ 7. ^ + Ci, (34) 

A- 2 2 A- 2 {A-2) ^B^^ +AC{A-2) ^ ' 

but we do not know how to obtain an "explicit" solution of this ode i.e. a solution of the form H = H{t). Possibly 
the most general solution to this equation may result unphysical as we have pointed out in the case of a perfect fluid 
(see the appendix of ([11])). 



1. Invariant solution 



In this case we can try to find a particular solution of eq. (23) through the invariant solution that induces the 
scaling symmetry X2 = [ai, —aH] . In such case we find that 



dt _ dH 



H=-, aeR, (35) 



which satisfies the eq. (23) iff a = a{uj), i.e. 



_ B±^B'^+ACiA~2) (3 + V3)±y^l2-9g^ 

2C |(^_l) + 3v/3(a; + l)' ^ ' 

In this case, it is observed that 

P=-H^ U = -^H'-^H'=(l-a)p, (37) 

therefore 



U»p=^U = >cp, x=y--ajeR^, (38) 

i.e. we have found that the viscous pressure and the energy density have the same order of magnitude and hence we 
can define a new equation of sate 11 = >cp where for physical reasons the numerical constants >c must be negative. 
Note that if w = 1, then >c = — |. At the same result R. A. Daishev and W. Zimdahl have arrived ( [X]) through a 
very different way. Nevertheless in our solution it is observed that >r can take other values. 

Note that both quantities have the same dimensional equation i.e. [11] = [p] and that for this reason under a scaling 
transformation the dimensionless quantity S must be remain constant (see the pioneering work in this field of D. M. 

Eardly [12], and the latter of K. Rosquits and R. Jantzen [13], and J. Wainwright [14]). Under the action of a similarity, 
each physical quantity (f> transforms according to it dimension q under scale transformations i.e. changes of the unit 
of length. Thus if unit of length L transforms as L — > XL then cj) — > A''0. This means that dimensionless quantities 
are invariant under a similarity transformation. Dimensionless quantities are therefore spacetimes constants. This 
implies that two quantities with the same dimensions, for example 11 and p or p and p are related through equations 
of state of the form 11 = >cp or p = top since the ratios S. or ^ must be constants. Furthermore, as J. Wainwright 
have pointed out, spacetimes admitting transitively self-similarity groups correspond exactly to the exact power law 
solutions as we have found. 

In the same way we can try to find a particular solution of eq.(28) that induces the symmetry X — xdx — 2ydy. 
Therefore we find that 



dx dy a 
X 2y ^ 



(39) 



is a solution of eq. (28) iff 



-B±JB^+4C{A-2) 
a = 



2C 

note that d = —a. Now tanking into account the change of variables (^y{x) = -ff7,x = iJ) it is found that 



where a is given by eq. (36). 



2. The case with 7 = 1/2, scale invariant solution and uj = 1, stiff matter. 
We shall study the equation 

H" - ^H-^ {H'f + (3 + Vs) HH' + 3V3i/3 = 0, 
which is a special case of eq. (23). (canonical variables) 



y = 3\/3.T3y^ + ( 3 + \/3 ) a;y2 



3y 

2a;' 



which is an Abel ode. The following change of variables brings us to obtain the next new ode 



s{r) = ~2hix,r = yx^ =^ x = e 



-l/s(r)' 



r (1 + 3^/3r2 + 2 (3 + VS) r) ' 

and which solution is: 

_ 2\/31n(6r + V3-I) 8\/31nr 2\/3 In (6r + \/3 + 3) 

V3-I (V3-l)(V3 + 3) V3 + 3 

2V31n(6y2;2 + x/3 - 1) 8^3 In (^/.t^) 2^3 In (Byx^ + x/3 + 3) 



-21nx 



V3-1 (\/3 - 1) (V3 + 3) \/3 + 3 

and hence 



-21nH 



2V31n(6f^ + V3- 1) 8%/31n(fi) 2^3 In (ef^ + V3 + 3 



V3-1 (V3-l)(V3 + 3) V3 + 3 



The invariant solution that we can find in this case is: 



which satisfies the eq. (42) iff 



t 



(3 + V3) ± V12 ^ f i + ^ = 0.7886751351, 
" 6v^ 1 -1 + ^ = 0.122008468. 
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IV. THE CASE n = xp. 

Since under a scale transformation we have found that the viscous parameter must be 7 = 1/2, and that in such 
case n — Kp^ now we are interested in studying the inverse way i.e. if under the hypothesis 11 ^ xp, the resulting 
differential equation remains scale invariant. For this purpose we rewrite the field eq. (3-6) under the assumption 
n = Kp. In this way the field equations may be expressed by the following ode: 

p"=P—-Al3p^p' + Bp\ (52) 
P 



with (1 — 7) = and 



■7 



^ ^ ^ K(a; + l + >.)(a; + l)(6 + 3>.) ^ ^ 

2SD 2x ' L J L J 

2x ( W 

w^{a^H), q; = (w + 1), — , ^=1 



6 + 3x V 2 / 2(tj + 1) 

The Lie analysis of equation (52) brings us to obtain the following system of pdes 

ipp + p-\p = Q. (53) 

_ ^^p-i + (^^^ _ 26p) + ^pAi^p^ = 0, (54) 

(277tp - 6t) + itAlip^ - -iipBp^ - ??t2p-i + 77^/32/"' = 0, (55) 

-^2Bp + r;„ + i^pEp" - 2^tBp' + ijtApp^' = 0, (56) 
we solve eqs. (53-56), finding that this system only admits the symmetry 

Xi=dt. (57) 
Now, if we try to check if the system admits a scaling symmetry 

X2 = afitdt + apdp, a e M, (58) 

we see with the help of eq. (56) that 

~2Aap^ + Aap^ + 2af3Ap'^ = 0, (59) 

finding in this way that 

-1 + 2/3 = 0^/3 = i, (60) 



where (3 = (1 — 7) , that is to say 7 = 5 



2 ■ 



^(p,t) = -|< + &, i^{p,t) = ap, (61) 

Therefore as we can see the ode only admits a single symmetry Xi = dt, but if we impose that the ode admits a 
scaling symmetry, we have seen that this is only possible if 7 = ^. Therefore iff 7 = ^ the ode admits two symmetries: 



Xi=du X2^tdt-2pdp, [Xi,X2]^Xi, (62) 

where X2 is the generator of the scaling group. Hence we can conclude that the assumption 11 = >cp does not imply 
that the resulting field equations must be scale invariant, this is only possible if 7 = i. 

The symmetry Xi brings us to obtain through the canonical variables the following Abel ode: 



where 



y' = -Bx^y^ + Af3x^y^ - -, (63) 

X 



y^-- (64) 
p 
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A. Equation (52) with 7 = 1/2. 

The equation (52) with 7 = 1/2 yields 

p"=P—-^^pp' + Bp\ (65) 
p z 

where [A^] = [B] . 

The symmetry Xi brings us to obtain through the canonical variables the following Abel ode: 

y' = -Bx^y^ + ^V^y^ (66) 

where 

We would like to point out that eq. (66) admits the following scaling symmetry 

3 

X = xdoo- -ydy, (68) 

which induces the following change of variables 

r = 2/x3/2, s(r)=lnx => y = ^51^. x = e'^''\ (69) 
in such a way that eq. (66) yields 



r (1 + rA - 2r2B) ' 

where the solution of eq. (70) is: 



2 arctan h 1 , 

s(r) = 21nr- In (l + rA-2r2B) + VVA-+ 8By ^^^^^ ^^^^ 



A-irB 



^/WTsB 

hence 



2 arctan h \ ^)^2+r" 



2 arctan h 



(70) 



Ina- = 2 In (yx^f'^ - In (l + (yx^f') A - 2 {y^x^) b) + ^ ^2 ^ ^ ^^^^ 

and tanking into account the change of variables (67) yields 



which is a quadrature, but unfortunately we do not know how to obtain an "explicit" solution for this ode as in the 
above case. 



1. Invariant solution 



The invariant solution is obtained for a 7^ For this value of a eq. (65) admits a single symmetry 

£. {t, P) = atdt, 11 {x, y) = -2apdp, 



(74) 
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the knowledge of one symmetry X might suggest the form of a particular solution as an invariant of the operator X 
i.e. the solution of 

this particular solution is known as an invariant solution (generalization of similarity solution) . In this case 



where 



2 / 1 AB + A^±A^{8B + A^) 



A = 2(. + l), 5^(^^±l±if)(^±l)(6±M, (77) 

2h 



with k.^ ~ K ~ \ and making w = 1 



, (2 + >^) 



2 



^ = 4, B = 3^ (78) 

A particular solution of eq. (66) may be found by taking into account the symmetry X — xdx — ^ydy In this case 

3— = -2^=^y=^, (79) 
and taking into account the change of variables {x = p^y = it is founded the already known solution p = pot~'^. 



B. The General solution. 



In this case, the assumption 11 = xp allows us to obtain a complete solution to the field equations (3-6). 
If we take into account eq. (6) with the assumption 11 = it yields: 

' [a + K) 2 {a + K) 2 

where H has been obtained from eq. (5) and follows the relationship 



H^-TT^ r-. (81) 



Simplifying eq. (80) it yields 



which trivial solution is: 



where 



a + X 2(a + >f) 2 J ^ """^ 

p^por^, (83) 



= -xk-'p''-^ (82) 



Po = 7 (84) 



<w 

a+K 2{a+K) 2 

In this way we have obtained a complete solution for the field equations (3-6) and valid for all value of 7. It is 
obvious that when 7 = 1/2 then we recover our previous solution p = pot^^- 
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V. A PEDESTRIAN METHOD. 



Ill this section we would like to show how dimensional Analysis works in order to obtain the same results but in 
a trivial way (see for example [15], [16], [17], and [18]). In the first place we will show how by writing the field 
equations in a dimensionless way we can determine the exact value of the parameter 7 which remains the equations 
scale invariant. In second place we would like to show how to solve some of the different odes that have arisen in this 
paper. 

Writing the field equations (3-6) in a dimensionless way and taken into account the following equations of state (8) 
it is obtained the following tt — monomia (see [19]): 



Gpt'' 



TTl = 



TTe = - = tH, 



7r2 = 



Ope 



TP 



it is observed that from Try 



and TTg 



we obtain 

TP 



n 

7r4 = — 
P 



T 



P 

DppP p 



(85) 
(86) 



(87) 



and from 713 and TTg 



Gpt-" 



k^p'^ 



1/(7-1) 



therefore 



c 

Gt 



1/(7-1) 



G 



j-b+2 



(89) 



where h — 1/(7— 1), it is observed that the only case and only for this, 7 — 1/2, wc obtain the relationship fc^ = (? jG. 
If 7 7^ 1/2 the "constants" G or c must vary or wc need to impose the condition Gjt? = const, (if both constants 
vary) if we want our equations to remain scale invariant as we have showed in an earlier paper (see [-()]) where we 
studied a viscous model with G time- varying. In such work we arrived to the conclusion that if 7 = 1/2, G must be 
constant in spite of considering it as a function that vary on time t, since we were only interested in the self-similar 
solution of that model. 

Now we go next to solve some of the differential equations that have arisen in this paper through the Dimensional 
technique. 

We begin studying eq. (24) i.e. 



H" - AH-^ {H' f + BHH' + CH^ = 



(90) 



which verifies the principle of dimensional homogeneity taking into account the dimensional base B = {T} . In this 
case we trivially arrive to the solution H oc since [H] = T~^. Note that D.A. (Pi theorem) does not understand 
numerical constants only of orders of magnitude. 
In second place we study eq. (28) 



y' = Cx^y^ + Bxy 



A 



y 



(91) 



with respect to the dimensional base B = {T} . This ode verifies the principle of dimensional homogeneity with respect 
to this dimensional base. Note that [y\ = [■^] = T^, and [x] = [H] = hence [y'\ = T^. Therefore rewriting the 
equation in a dimensionless way we find that y oc x"^ 

But if we study this equation with respect to the dimensional base B = {X, Y} , we need to introduce new 
dimensional constants that make that the equation verifies the principle of dimensional homogeneity 



y' = aCx^y^ + [3 Bxy 



A- 



(92) 
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where [a"'^^^] 



hence 





y 


/3 


X 


X 







1 


Y 


1 


-1 










(93) 



As we can see we have obtained the same solution than in the case of the invariant solution. This is because the 
invariant solution that induces a scaling symmetry is the same as the obtained one through the Pi theorem. 

We would like to emphasize that D.A. brings us to obtain change of variables (c.v.) (see [21] for more details) which 
allows us to obtain odes simplest than the original one. in this case, it is observed that [f3] = X~'^Y~^ in such a way 
that we have the c.v. 



X ^ t. 



therefore eq. (92) yields: 



{t = X, u{t) = Px'^y) =^ 

tu =■ u {y? + u + l) , 



y 



and hence 



1/9 s •\/3 

hit ^ — In [W + ii + 1 j H arctan 

2 3 



Inu + Ci = 0, 



in the original variables it yields 

Xnx+\ln({ax^yf + ax^V+l) + 
In the same way we can study eq. (65) . 



V3 



arctan 



P ^ 



where = [B] = LAI ^, with respect to the dimensional base B = {L, M, T} . Therefore it is found 





P 


B 


t 


L 


-1 


1 





M 


1 


-1 





T 


-2 





1 



p oc 



Bt^' 



To end we study eq. (66) 



y'^^BxV + ^V^y'-^, 

2 x 



where [^^] = [B] — X ^, with respect to the dimensional base B ~ {X,Y} . Therefore it is found 





y 


B 


X 


X 





-2 


1 


Y 


1 


-1 






y c< 



as already we know. 



(94) 
(95) 
(96) 



( Q (ax'y) + 1)^)- 111 i^^^y) +Ci=0. (97) 



(98) 



(99) 



(100) 



VI. CONCLUSIONS 



In this paper we have studied the possible symmetries that admits a flat FRW model filled with a bulk viscous 
fluid. Using the Lie group method we have tried to find an adequate equation of state for the viscous parameter as 
well as for the viscous pressure. Therefore we conclude that the field equations remain scale invariant iff 7 = 1/2 
and that for this parameter it is found that 11 = >cp. But that the hypothesis 11 ~ Kp does not imply that the field 
equations remain scale invariant, this only occurs if 7 = 1/2. Furthermore, the assumption 11 = yep brings us to obtain 
a complete solution to the field equations valid for all 7. 
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